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(N 1 INTRODUCTION 

• • The measurement of the Rotation Curves (RCs) of disk 
. — . galaxies is a powerful tool to investigate the nature of Dark 
Matter (DM), including its content relative to the bary- 
» I onic components and their distributions. In particular, dwarf 
C3 galaxies are good candidates to reach this aim as their kine- 
matics are generally dominated by the dark component, 
down to small galactocentric radii (Persic et al. 1996; Gen- 
tile et al. 2004; de Blok et al. 2001; Oh et al. 2008; Salucci et 
al. 2008). This leads to a reliable measurement of the dynam- 
ical contribution of the DM to the RC and hence of its den- 
sity profile. Therefore, a dwarf galaxy like the Orion dwarf 
provides us with an important test as to whether DM density 
profiles arising in A Cold Dark Matter (ACDM) numerical 
simulations (Navarro et al. 1996) are compatible with those 
detected in actual DM halos around galaxies. Let us com- 
ment that NFW profile arises from pure N-body DM simu- 
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ABSTRACT 

Dwarf galaxies are good candidates to investigate the nature of Dark Matter, because 
their kinematics are dominated by this component down to small galactocentric radii. 
We present here the results of detailed kinematic analysis and mass modelling of the 
Orion dwarf galaxy, for which we derive a high quality and high resolution rotation 
curve that contains negligible non-circular motions and we correct it for the asymmet- 
ric drift. Moreover, we leverage the proximity (D — 5.4 kpc) and convenient inclination 
(47°) to produce reliable mass models of this system. We find that the Universal Rota- 
tion Curve mass model (Freeman disk + Burkert halo -I- gas disk) fits the observational 
data accurately. In contrast, the NFW halo + Freeman disk + gas disk mass model 
is unable to reproduce the observed Rotation Curve, a common outcome in dwarf 
galaxies. Finally, we attempt to fit the data with a MOdified Newtonian Dynamics 
(MOND) prescription. With the present data and with the present assumptions on 
distance, stellar mass, constant inclination and reliability of the gaseous mass, the 
MOND "amplification" of the baryonic component appears to be too small to mimic 
the required "dark component" . The Orion dwarf reveals a cored DM density distribu- 
tion and a possible tension between observations and the canonical MOND formalism. 

Key words: Dark Matter; Galaxy: Orion dwarf; Mass profiles 



lations. It is well known that, as effect of the baryonic infall 
in the cosmological DM halos and of the subsequest process 
of stellar disk formation, shallower profiles of the DM halo 
may arise (see Governato et al. (2010, 2012); Maccio' et al. 
(2012)). 

Recent studies of the RCs of dwarf galaxies have tested 
the NFW scenario. It is now clear that kinematic data are 
better fitted by a DM halo with a constant density core (e.g. 
Borriello & Salucci (2001); van den Bosch & Swaters (2001); 
de Blok et al. (2001); Weldrake et al. (2003); Gentile et al. 
(2004); Simon et al. (2005); Oh et al. (2008)), than by one 
that is centrally peaked. One specific example is DDO 47, 
whose velocity field is clearly best fitted if the DM halo is 
cored; moreover, its (small) detected non-circular motions 
cannot account for the discrepancy between data and the 
NFW predictions (Gentile et al. 2005). 

The present investigation examines the DM content of the 
Orion dwarf galaxy. This nearby system harbors an extended 
H I disk, and thus provides us with an important test of the 
above paradigm. As we show below, the Orion dwarf is one 
of the few known galaxies whose kinematics unambiguously 
point towards a cored profile. This system is thus critically 
important for investigating the nature of the DM particle 
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Figure 1. Comparison of the Hi column density distribution (a), the intensity- weighted velocity field (b), and the intensity- weighted 
velocity dispersion (c); the beam size is 20". The contour in (a) is at the 10 21 cm -2 level; the contours in (b) show velocities between 
320 and 420 km s — 1 , separated by 10 km s — , and the thick white line corresponds to the 400 km s _1 contour; the contours in (c) are 
at the 5, 10, 15 km s _1 levels. 



and of the evolution of DM halos. 

MOND accounts for the evidence that RCs of spiral galaxies 
are inconsistent with the corresponding distribution of the 
luminous matter (Milgrom 1983a,b). Rather than postulat- 
ing the existence of a dark halo made by massive collision- 
less elementary particles, this scenario advocates that the 
gravitational force at low accelerations leaves the standard 
Newtonian regime to enter a very different one. Historically 
MOND has generally been successful in reproducing the RCs 
of spiral galaxies with only the (observed) luminous matter 
(e.g. Sanders (1996); Sanders & Verheijen (1998); Sanders 
& Noordermeer (2007)). However, cases of tension between 
data and the MOND formalism do exist (Gentile et al. 2004). 
It is important to stress that in order to derive the DM den- 
sity profile or to test the MOND formalism, we must know 
the distribution of the ordinary baryonic components, as 
well as have reliable measurements of the gas kinematics. 
For the Orion dwarf, 21-cm Hi surface brightness and kine- 
matics have recently been published (Cannon et al. 2010): 
their analysis provides a high quality, high resolution RC, 
that, in addition, can be easily corrected for asymmetric 
drift and tested for non-circular motions. This galaxy is a 
very useful laboratory in that a simple inspection of the 
RC ensures us that it shows a large mass discrepancy at 
all radii. Moreover, the baryonic components are efficiently 
modeled (i.e., no stellar bulge is evident and the stellar disk 
shows a well-behaved exponential profile, see Vaduvescu et 
al. (2005)). The distance to the galaxy, which is critical for 
an unambiguous test of MOND (Sanders & McGaugh 2002), 
is estimated to be 5.4±1.0 Mpc (Vaduvescu et al. 2005). It 
is important to stress that the distance of the Orion dwarf 
remains a significant source of uncertainty. Vaduvescu et 



al. (2005) estimate the distance using the brightest stars 
method. The intrinsic uncertainty in this technique may al- 
low a distance ambiguity much larger than the formal errors 
estimated by Vaduvescu et al. (2005), because in their work 
this method yields a scatter as large as 50% in distance. 
Finally, the system's inclination (47°) is kinematically mea- 
sured (see Section (3.1)) and is high enough to not affect the 
estimate of the circular velocity. The properties described 
above make the Orion dwarf galaxy an attractive candidate 
to determine the underlying gravitational potential of the 
galaxy. 

This paper is organized as follows. In Sec. 2 we present 
the stellar surface photometry. In Sec. 3, the Hi surface 
density and kinematics data are presented and discussed; we 
also provide the analysis of possible non-circular motions of 
the neutral gas. In Sec. 4 we model the RC in the stellar 
disk using a cored/cusped halo framework. In Sec. 5 we test 
the Orion kinematics against the MOND formalism. Our 
conclusions are given in Sec. 6. 



2 STELLAR PHOTOMETRY 

Following the discussion in Cannon et al. (2010), the under- 
lying stellar mass in the Orion dwarf is estimated using the 
near-infrared (IR) photometry (J and Ks bands) presented 
by Vaduvescu et al. (2005). Those authors find (J — Ks) = 
+0.80 and a total Ks magnitude of +10.90. When compar- 
ing to models (see below) we assume that the color differ- 
ence between K and Ks is negligible; further, we assume 
Lk,© = +3.33 (Cox 2000; Bessel 1979). Accounting for ex- 
tinction, the total K-band luminosity of the Orion dwarf 
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is ~3.5 x 10 8 L®. The mass of the stellar component was 
estimated by Cannon et al. (2010) to be (3.7±1.5) x 10 8 
M®. The stellar surface brightness profile is well fitted by 
an exponential thin disk, with a scale length of Rd= 25" ± 
1" (equivalent to 1.33 ± 0.05 kpc at the adopted distance). 
Moreover, there are no departures from an exponential pro- 
file that would be indicative of a prominent central bulge. 



3 HI SURFACE DENSITY AND KINEMATICS 

H I spectral line imaging was acquired with the Very Large 
Array and presented in Cannon et al. (2010). We refer the 
reader to that work for a full discussion of the data handling, 
and we summarize salient details here. The final data cubes 
have a circular beam size of 20", with a 3 a Hi column 
density sensitivity of Nhi= 1.5 xl0 19 cm~ 2 . The first three 
moment maps (i.e. the integrated Hi intensity, the velocity 
field, and the velocity dispersion) are shown in Figure (1). 
The neutral gas disk of the Orion dwarf shows rich morpho- 
logical and kinematic structure at this physical resolution. 
The outer disk contains tenuous Hi gas, but column densi- 
ties rise above the 5 x 10 20 cm -2 level at intermediate radii. 
There is plentiful high-column density (>10 21 cm" 2 ) Hi 
throughout the disk. The more or less parallel iso-velocity 
contours at inner radii are indicative of linear rotation (al- 
though almost certainly not solid body) and the curving of 
the outer contours suggests that the outer rotation curve has 
a fairly constant velocity. The outer disk contours show no 
evidence for a decrease in rotational velocity at large radii. 
In the central regions of the disk, however, some H I "holes" 
or "depressions" manifest a pronounced kink in these con- 
tours (consider the contours at 370 ± 20 km s _1 ). The inten- 
sity weighted velocity dispersion averages to ~7-8 km s _1 
throughout the disk, although the innermost regions show 
dispersions above 10 km s _1 . 

The total Hi flux integral, proportional to the Hi disk 
mass, was found to be 50.3 ±5.1 Jykm s _1 , a value some- 
what lower than the single-dish flux measure of 80.6±7.72 
Jykm s _1 by Springob et al. (2005); the difference may arise 
from the lack of short interferometric spacings that provide 
sensitivity to diffuse structure. The total H I mass is found 
to be Mm = (3.5±0.5) x 10 8 M Q . After applying the usual 
35% correction for Helium and molecular material, we adopt 
M gas = (4.7±0.7) x 1O 8 M as the total gas mass. 
In Figure (2) we plot the 10" /20" resolution H I surface den- 
sity, throughout the gas disk. A simple fit (valid out to the 
last measured point and for the scope of this work) yields: 

-0.263r 3 + 1.195r 2 + 3.094r+ 18.549 „„ , 2 

^ Hl(r) = 0.154r3 - 1.437r + 6.703 M ° /pC ' 

(1) 

where r is in kpc. The related fitting uncertainty on j-Lm{r) 
is about 20%. Figure (2) shows that the Hi surface density 
rises from the center of the galaxy, reaches a maximum, and 
then declines exponentially. At the last measured point, i.e. 
out to ~ 7 kpc, the profile has almost (though not com- 
pletely) reached the edge of the Hi disk and rapidly con- 
verges to zero. Note that, in Newtonian gravity, the outer 
gaseous disk contributes in a negligible way to the galaxy 
total gravitational potential. 



3.1 The Circular Velocity 

The channel maps of the Orion dwarf provide evidence of 
well-ordered rotation throughout the H I disk (see Cannon 
et al. 2010). The intensity- weighted-mean velocity field (Fig- 
ure (lb)) exhibits symmetric structure in the outer disk. 
Twisted iso-velocity contours at inner radii coincide with 
the H I holes near the centre of the disk. The disk is there- 
fore dominated by circular motion. The RC of the galaxy 
was derived by fitting a tilted ring model to the intensity- 
weighted-mean velocity field using the GIPSY task ROTCUR. 
The routine carries out a least-squares fit to Vi os , the line 
of sight-velocity. To derive the best-fitting model, an iter- 
ative approach was adopted in which the various combina- 
tions of the parameters were fitted. The final RC was ex- 
tract by fixing all other parameters. The receding and ap- 
proaching sides of the galaxy were fitted separately. The 
best fitting parameters are i= (47 ± 3)°, P. A. = (20 ± 2)°, 
V sys = 368.5 ± 1.0 km/s, and (02000,^2000) = (05:45:01.66, 
05:03:55.2) for the dynamical centre. We have realized that 
the inclination is not dependent on the radius, and the fit is 
shown in Figure (3). Its weighted value is (46.8±0.14)°. No- 
tice that because the errors reported by GIPSY/ROTCUR 
include only errors on the fits and systematics are not in- 
cluded, the 3° error estimate comes from attempting the 
ROTCUR fits in various orders (e.g., holding each variable 
fixed in turn). 

The resulting RC is shown in Figure (4). Notice that in this 
object the disk inclination is determined kinematically and 
therefore it is quite accurate. No result of this paper changes 
by adopting different values of i, inside the quoted errorbar. 

The second-order moment map for the galaxy is shown 
in Figure (lc). Throughout most of the disk, the velocity 
dispersion is roughly constant at a ~ 7±2 km/s, with a 
more complex behaviour near the galaxy centre and at the 
outermost radii. This velocity dispersion estimate allows us 
to derive the asymmetric drift correction to the RC yielded 
by the tilted ring model. The observed rotation velocity, 
Vrot, is related to the circular velocity V c via 



V c \r) 



V(r) 



dlog /im(r) dlog cr 2 (r) 



dlog r 



dlog r 



(2) 



From an examination of Figure (4) it is clear that the 
Vrot and V c profiles differ by less than 1%. Throughout this 
paper, we use the latter for the purposes of mass modelling. 
We notice that in very small dwarfs this correction is not 
negligible (V ro t ~ <r) and it introduces an uncertainty in the 
analysis, e.g. Begum & Chengalur (2004). 

In summary, the Orion dwarf RC has a spatial resolu- 
tion of 0.26 kpc (i.e. 0.2 Rd), and extends out to 5.1 Rr>. 
The uncertainties on the RC are few km/s and the error on 
the RC slope dlogV/dlogR < 0.1. 

Is the circular velocity given by eq. (2) a proper estimate 
of the gravitational field? To further investigate the presence 
of non-circular motions within the Hi disk that jeopardize 
the kinematics, we carried out a harmonic decomposition of 
the intensity-weighted velocity field to search for any signifi- 
cant non-circular components. This test is necessary in that 
the undetected presence of non-circular motions can lead to 
incorrect parametrization of the total mass distribution. 
Following Schoenmakers et al. (1997), the line-of-sight ve- 
locities from the H I velocity field are decomposed into har- 
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Figure 2. Radially averaged Hi mass surface density profiles of the Orion dwarf, created by averaging Hi emission in concentric rings 
emanating from the dynamical center found in our RC analysis. The solid/dotted lines were created from the 10" /20" resolution images. 
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Figure 3. Inclination best fit (thick Black line) for high resolution (10") data (filled Blue points). Errorbars are plotted and in most 
cases are smaller than the symbol size. The errors reported (by GIPSY/ROTCUR) include only errors on the fits, while the systematic 
error considered (3°) is shown as a Green errorbar. 



monic components up to order N = 3 according to 

3 

Vios = Vsys + ^2 c m cosmfl + s m sinra#, (3) 

m — 1 

where Vsys is the systemic velocity, c m and s m are the mag- 
nitudes of the harmonic components, m the harmonic num- 
ber, and the azimuthal angle in the plane of the galaxy. 
The GIPSY task RESWRI was used to carry out the decom- 
position by fitting a purely circular model to the velocity 
field, subtracting it from the data, and then determining 



from the residual the magnitudes of the non-circular com- 
ponents. The tilted ring model fitted by RESWRI had its 
kinematic centre fixed to that of the purely circular tilted 
ring model used to derive the RC above. The position angles 
and inclinations were fixed to constant values of 20 and 47° , 
respectively. 

The parameters of the best-fitting model are shown in 
Figure (5). Adjacent points are separated by a beam width 
in order to ensure that they are largely independent of one 
another. We argue that because the standard tilted ring 
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Figure 4. Rotation velocity corrected for asymmetric drift (filled black points), raw velocity data with error bars (filled red boxes) and 
the asymmetric drift correction (filled blue triangles). 



model has fewer free parameters than the model incorpo- 
rating the higher order Fourier components, it is not as es- 
sential to space the points on the rotation curve by a full 
beam width. Then in this model only 16 points are consid- 
ered instead of the 26 points used in fitting the RC. 
At inner radii the inferred non-circular motions are not neg- 
ligible, but this is almost certainly due to the fact that the 
Hi distribution over this portion of the disk is irregular, 
being dominated by the large central Hi under-densities. 
The harmonic components of the outer disk are, instead, 
reliable and demonstrate the gas flow to be dominated by 
circular kinematics. The circular velocity so obtained well 
matches that found by means of the tilted ring model pre- 
sented above. The amplitudes of C2 and S2 are too small to 
hide a cusp inside an apparently solid body RC (as suggested 
by Hayashi & Navarro (2006)). These results provide further 
decisive support for the use of V ro t of the Orion dwarf as a 
tracer of its mass distribution. 



4 MASS MODELING 

We model the Orion dwarf as consisting of two "luminous" 
components, namely the stellar and the gaseous disks, em- 
bedded in a dark halo. The stellar component is modelled as 
an exponential thin disk (Freeman 1970) with a scale length 
of 1.33 kpc. Any bulge component is assumed to be negli- 
gible in terms of mass. The dynamical contribution of the 
gas to the observed RC is derived from the H I total inten- 
sity map. A scaling factor of 1.33 is incorporated to account 
for the presence of Helium and other elements. For the dark 
halo we consider two different parametrizations of the mass 
distribution: an NFW profile (Navarro et al. 1996) and the 
cored profile of the Halo Universal Rotation Curve (URCH) 
(Salucci et al. 2007). It is well known that the NFW profile 
is one outcome in numerical simulations of cold dark matter 
structure formation, whereas the cored profile (an empirical 
result), by design, fits the broad range of RC shapes of spiral 
galaxies. 



4.1 Mass Models: stellar disk + dark matter halo 

The RC is modelled as the quadrature sum of the RCs of 
the individual mass components: 



+ V 



(4) 



For the cored halo parametrization we adopt the URCH pro- 
file: 



6.4 



Po4 



In 1 + 



'0 



arctan 



r 



1 , / - r 

-In 1 + ^ 



(5) 



where the disk mass, the core radius ro and the central halo 
density po are free parameters. 

It is evident that this model yields a total RC that fits 
the data extremely well (see Figure (6) left panel), with 
best-fitting parameters of ro = (3.14 ± 0.32) kpc, M_d = 
(3.5 ± 1.8) x 10 s M and p = (4.1 ± 0.5) x 10~ 24 g/cm 3 . 
More accurate statistics is not necessary; the mass model 
predicts all the V(r) data points within their observational 
uncertainty. 

Notice that the derived value of the disk mass agrees 
with the photometric estimate discussed above. The cor- 
responding virial mass and radius of the DM halo are 
M vir = (5.2 ± 0.5) x 10 10 M (see eq. 10 in Salucci et al. 
(2007)) and R^ r = 95. 5±| kpc (Eke et al. 1996), respec- 
tively. We note that the Orion dwarf has a mass 20 times 
smaller than that of the Milky Way, with the DM halo domi- 
nating the gravitational potential at all galactocentric radii. 
The baryonic fraction is fb = (Md + M sas )/M V ir = 0.016, 
while the gas fraction is M gas /M„ 4r = 9 x 10 -3 . 

The RC for the NFW dark matter profile is 



V NFW (r) = V v . 



9(c) 
'xg(cx) 



(6) 



where x — r/R vir , g(c) — [ln(l + c) — c/(l + c)] and 
c(M vir ) = 9.60 (Ai^/lO^/i^M©)" 0,075 is the concentra- 
tion parameter (see Klypin et al. (2011)). We fitted the RC 
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Figure 5. Non-circular velocity components as derived by carrying out a harmonic decomposition of the intensity-weighted-mean velocity 
field. The c\ profile corresponds to the circular RC. The RC derived by fitting tilted ring model to the intensity-weighted-mean velocity 
field is shown as a solid blue curve. 




Figure 6. Left: URC model of Orion dwarf. The circular velocity (filled circles with error bars) is modeled (thick Red line) by a Freeman 
disk (Magenta line), a URCH halo (Green line), and the Hi circular velocity (Black line). Right: NFW model of Orion dwarf. The 
circular velocity data (filled circles with error bars) is modeled (thick Blue line) including a Freeman disk profile (Magenta line), a NFW 
halo profile (Orange line) and the Hi circular velocity (Black line). In both cases the values of the free parameters are reported in the 
text. 



of the Orion dwarf by adjusting M„,: r and Mu. The result- 
ing best-fit values are M vir = (2.5 ± 0.5) x 10 11 M Q and 
M D = (6.9 ±1.7) x 10 s M Q , but since x^ e d ~ 3 - 3 > Le - the fit 
is unsuccessful, the best-fit values of the free parameters and 
those of their fitting uncertainties do not have a clear physi- 
cal meaning. We plot the results in the right panel of Figure 
(6). The NFW model, at galactocentric radii r<2 kpc, over- 
estimates the observed circular velocity (see Figure (6) right 
panel) . 



5 ORION KINEMATICS IN THE MOND 
CONTEXT 

An alternative to Newtonian gravity was proposed by Mil- 
grom (1983a) to explain the phenomenon of mass discrep- 
ancy in galaxies. It was suggested that the true acceleration 
a of a test particle, at low accelerations, is different from the 
standard Newtonian acceleration, ajv: 



H{a/a ) 

where (i{a/ao) is an interpolation function and ao = 1.35 x 
10 _8 cm s -2 is the critical acceleration at which the transi- 
tion occurs (see Famaey et al. (2007)). For this we adopt the 
following form of the interpolation function (see Famaey & 
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Binney (2005)): 

1 + a/oo 

In this framework the circular velocity profile can be 
expressed as a function of a and of the standard Newtonian 
contribution of the baryons to the RC, Vi, ar =(V| 3 +Vg as ) 1 ^ 2 , 
obtaining for it 

V 2 - V 2 M ( "J 1 + V t^ + 1 \ (q) 

VMOND — Vbar\ r ) I ^ I , l»J 

(see Richtler et al. (2008)). Eq. (9) shows that in the MOND 
framework the resulting RC is similar to the no-DM stan- 
dard Newtonian one, with an additional term that works 
to mimic and substitute for the DM component (Sanders & 
McGaugh 2002). No result of this paper changes by adopt- 
ing the "standard" MOND interpolation function (see Zhao 
& Famaey (2006)). 

The best-fitting MOND mass model is shown in Fig- 
ure (7). The model total RC (cyan line) completely fails 
to match the observations. We fix the stellar mass Md at 
Md = 2.6 x 1O 8 M . If we let the disk mass becomes higher, 
covering the mass range estimated in Cannon et al. (2010), 
the fit is not even able to reproduce the RC at inner radii. 
Note that in the MOND formalism, the distance of the 
galaxy and the amount of gaseous mass are both crucial in 
deriving the model RC. To quantify the discrepancy of these 
observations with the MOND formalism, note that only if 
the Orion dwarf were 1.9 times more distant than the cur- 
rent estimate we would obtain a satisfactory fit to the RC 
(see Figure (8)). 



6 CONCLUSIONS 

The Orion dwarf galaxy is representative of a population of 
dwarfs with a steep inner RC that gently flattens at the edge 
of the gas disk. The observed kinematics imply the presence 
of large amounts of DM also in the central regions. We have 
used new H I observations of the Orion dwarf to analize its 
kinematics and derive the mass model. The derived RC is 
very steep and it is dominated by DM at nearly all galac- 
tocentric radii. Baryons are unable to account for the ob- 
served kinematics and are only a minor mass component at 
all galactic radii. 

We have used various mass modeling approaches in this 
work. Using the NFW halo, we find that this model fails to 
match the observed kinematics (as occurs in other similar 
dwarfs). We show that non-circular motions cannot resolve 
this discrepancy. Then we modelled the galaxy by assuming 
the URCH parametrization of the DM halo. We found that 
this cored distribution fits very well the observed kinemat- 
ics. Orion is a typical dwarf showing a cored profile of the 
DM density and the well-known inability of DM halo cuspy 
profiles to reproduce the observed kinematics. Finally, we 
find that the MOND model is discrepant with the data if 
we adopt the literature galaxy distance and gas mass. The 
kinematic data can be reproduced in the MOND formalism 
if we allow for significant adjustments of the distance and/or 
value of the gas mass. Let us point out that the present in- 
terferometric observations may miss some of the objects' Hi 



flux, although this may be limited in that the cubes do not 
have significant negative bowls. Obviously, for bigger values 
of the Hi mass, the distance at which the baryon compo- 
nents would well fit the data will also somewhat decrease. 
It is worth stressing that there is a galaxy distance (albeit 
presently not-favoured) for which MOND would strike an 
extraordinary success in reproducing the observed kinemat- 
ics of the Orion dwarf. 

The Orion dwarf has a favorable inclination, very reg- 
ular gas kinematics, a small asymmetric drift correction, a 
well-understood baryonic matter distribution, and a large 
discrepancy between luminous and dynamical mass. All of 
these characteristics make this system a decisive benchmark 
for the MOND formalism and a promising target for further 
detailed studies. Of particular value would be a direct mea- 
surement of the distance (for example, infrared observations 
with the Hubble Space Telescope would allow a direct dis- 
tance measurement via the magnitude of the tip of the red 
giant branch). 
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